A New Limit on Time-Reversal Violation in Beta Decay 
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We report the results of an improved determination of the triple correlation DP • (pe x pv) that 
can be used to limit possible time-reversal invariance in the beta decay of polarized neutrons and 
constrain extensions to the Standard Model. Our result is D = (—0.96 ± 1.89{stat) ± l.Ol(sys)) x 
10^*. The corresponding phase between gA and gv is 4>av ~ 180.013° ± 0.028° (68% confidence 
level). This result represents the most sensitive measurement of D in nuclear beta decay. 

PACS numbers: 24.80. ly, ll.30.Er, 12.15.Ji, 13.30.Ce 



The existence of charge-parity (CP) symmetry viola- 
tion in nature is particularly important in that it is nec- 
essary to explain the preponderance of matter over an- 
timatter in the universe '1]. Thus far, CP violation has 
been observed only in the K and B meson systems 0- 
0] and can be entirely accounted for by a phase in the 
Cabbibo-Kobayashi-Maskawa matrix in the electroweak 
Lagrangian. This phase is insufhcient to account for the 
known baryon asymmetry in the context of Big Bang 
cosmology [5|, so there is good reason to search for CP 
violation in other systems. As CP and time-reversal (T) 
violation can be related to each other through the CPT 
theorem, experimental limits on electric dipole moments 
and T-odd observables in nuclear beta decay place strict 
constraints on some, but not all, possible sources of new 
CP violation. 

The decay probability distribution for neutron beta de- 
cay, dW, can be written in terms of the beam polariza- 
tion P and the momenta (energies) of the electron Pe 
{Ee) and antineutrino Pi, {Ev) as 
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A contribution of the parity-even triple correlation DP • 
(Pe X p^) above the level of calculable final-state inter- 
actions (FSI) implies T-violation. The PDG average of 
recent measurements is D — (—4 ±6) x 10"'' 043 , while 
the FSI for the neutron are - 10"^ [13, Ell. Comple- 
mentary limits can be set on other T-violating correla- 
tions, and recently a limit on R has been published [l^ . 
Various theoretical models that extend the SM, such as 
left-right symmetric theories, leptoquarks, and certain 
exotic fermions could give rise to observable effects that 



are as large as the present experimental limits |13i]. Cal- 
culations performed within the Minimal Supersymmetric 
Model, however, predict D < 10""^ 



In the neutron rest frame, the triple correlation can be 
expressed as DP ■ (pp x Pe), where Pp is the proton mo- 
mentum. Thus one can extract D from the spin depen- 
dence of proton-electron coincidences in the decay of cold 
polarized neutrons. Our measurement was carried out 
at the National Institute of Standards and Technology 
Center for Neutron Research [l5[. The detector, shown 
schematically in Fig. [1] consisted of an octagonal array of 
four electron-detection planes and four proton-detection 
planes concentric with a longitudinally polarized beam. 
The beam, with a neutron capture fiuence rate at the 
detector of 1.7 x 10^ cm~^ s~*, was defined using a 
series of ^LiF apertures and polarized to > 91% (95% 
C.L.) by a double-sided bender-type supermirror [15'|. A 
560 f/T guide field maintained the polarization direction 
throughout the fiducial volume and a current-sheet spin- 
flipper was used to reverse the neutron spin direction 
every 10 s. The symmetric octagonal geometry was cho- 
sen to maximize sensitivity to D while approximately 
canceling systematic effects stemming from detector ef- 
ficiency variations or coupling to the spin correlations A 
and B [l^ . Each of the four proton segments consisted 
of a 2 x 8 array of silicon surface-barrier detectors (SBDs) 
with an active layer SOOmm^x 300 /xm. Each SBD was 
contained within an acceleration and focusing cell con- 
sisting of a 94% transmitting grounded wire-mesh box 
through which the recoil protons entered. Each SBD, 
situated within a field-shaping cylindrical tube, was held 
at a fixed voltage in the range — 25kV to — 32kV. The 
sensitive regions of the beta detectors were plastic scin- 
tillator measuring 50 cm by 8.4 cm by 0.64 cm thick, with 
photomultiplier tube (PMT) readout at both ends. This 
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FIG. 1: A schematic of the emiT detector illustrating the alternating electron and proton detector segments. The darker shaded 
proton detectors indicate the the paired-ring at z = ±10 cm. The cross section view illustrates, in an exaggerated manner, the 
effect of the magnetic field on the particle trajectories and average opening angle. A P2E3 coincidence event is shown. 
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FIG. 2: Intensity log plot of SBD-scintillator coincidence data 
showing proton energy vs delay time. Events near At — 
are prompt coincidences due primarily to beam-related back- 
grounds. 



thickness is sufficient to stop electrons at the decay end- 
point energy of 782 keV. The proton and beta detectors 
were periodicahy calibrated in situ with gamma and beta 
sources respectively. Details of the apparatus are pre- 
sented elsewhere [sl. TBI 17 . 

Data were acquired in a series of runs from October 
2002 through December 2003. Typical count rates were 
3s~^ and 100 s^^ for single proton and beta detectors, 
respectively, while the coincidence rate for the entire ar- 
ray was typically 25 s""'^. Of the raw events, 12% were 
eliminated by filtering on various operational parameters 
(e.g. coil currents) and by requiring equal counting time 
in each spin-fiip state. A beta-energy software threshold 
of 90 keV eliminated detection efficiency drifts due to 
changes in PMT gain coupled with the hardware thresh- 
old. This was the largest single cut, eliminating 14% of 
the raw events. A requirement that a single beta be de- 
tected in coincidence with each proton eliminated 7% of 
events. All cuts were varied to test for systematic effects. 



The remaining coincidence events were divided into 
two timing windows: a preprompt window from -12.3 fis 
to -0.75 /iS that was used to determine the background 



from random coincidences, and the decay window from 
-0.5 /is to 6.0 /iS as shown in Fig. [5] The recoil proton 
has an endpoint of 750 eV. On average it is delayed by 
^ 0.5 ^s. The average signal-to-background ratio was 
~ 30/1. The energy- loss spectrum produced by mini- 
mum ionizing particles in 300 /im of silicon is peaked at 
approximately 100 keV and, being well separated from 
the proton energy spectrum, yielded an estimated con- 
tamination below 0.1%. The final data set consisted of 
approximately 300 million accepted coincidence events. 

A detailed Monte Carlo simulation was used to esti- 
mate a number of systematic effects. The program Pene- 
lope J which has been tested against data in a variety 
of circumstances of relevance to neutron decay jl^ j was 
embedded within a custom tracking code. All surfaces 
visible to decay particles were included. The Monte Carlo 
was based on the measured beam distribution upstream 
and downstream of the fiducial volume 15| and incorpo- 
rated the magnetic field and electron energy threshold. A 
separate Monte Carlo based on the package SIMION , 
incorporating the detailed geometry of the proton cells, 
was used to model the proton detection response func- 
tion. 

Achieving the desired sensitivity to D in the presence 
of the much larger spin- asymmetries due to A and B 
depends critically on the measurement symmetry. To the 
extent that this symmetry is broken, corrections must be 
applied to the measured result. These corrections are 
listed in Table Hand are discussed below. To extract D, 
coincident events were first combined into approximately 
efficiency-independent asymmetries 



(2) 



where N^'^' is the integrated number of coincident events 
in proton detector i = 1...64, beta detector j — 1...4, 
with neutron spin + (— ) aligned (anti-aligned) with the 
guide field. For uniform polarization, P, the asymme- 
tries, lyP'^^j , can be written in terms of decay correlations 
as 
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wP'^^- « P . {A-k^f + B'K^b' + ijkS'' ), (3) 

where the K's are obtained from Eqn. [1] by integrat- 
ing the normahzcd kinematic terms over the phase space 
of the decay, the neutron beam volume, and the accep- 
tance of the indicated detectors ^.8|. Ka oc (pe/Ee) and 
Kb oc [Tpy/Eu) are primarily transverse to the detector 
axis but have roughly equal longitudinal components for 
coincidence events involving the two beta detectors op- 
posite from the indicated proton detector (E2 and E3 for 
P2 as shown in Fig.[T]). The K£)'s, however, are primarily 
along the detector axis and are opposite in sign for the 
two beta detectors. Thus for each proton detector we 
can choose an appropriate combination of detector pairs 
that is sensitive to the D-correlation but that largely can- 
cels the parity-violating A and B correlations. One such 
combination is 



TABLE I: Systematic corrections and combined standard un- 
certainties (68% confidence level). Values should be multi- 
plied by 10"". 



Source 


Correction 


Uncert ainty 


Background asymmetry 


Qa 


0.30 


Background subtraction 


0.03 


0.003 


Electron backscattering 


0.11 


0.03 


Proton backscattering 





0.03 


Beta tiireshold 




n 1 n 




—0.29 


0.41 


Beam expansion, magnetic field 


-1.50 


0.40 


Polarization non-uniformity 




0.10 


ATP - misalignment 


-0.07 


0.72 


ATP - Twist 


Qa 


0.24 


Spin-correlated fiux'' 


Qa 


3 X 10-*^ 


Spin-correlated pol. 


Qa 


5 X 10"" 


Polarization*^ 




0.04^* 


Kd" 




0.05 


Total systematic corrections 


-1.68 


1.01 



yP. ^ _(^yjP.en _^P.e^^^ (4) 

where cr and cl label the electron-detector at approxi- 
mately 135° giving a positive and negative cross-product 
Pp X pe respectively (P2E3 vs P2E2 as shown in Fig.[T|). 
Proton cells at the detector ends accept decays with 
larger longitudinal components of K^i and are more sen- 
sitive to a range of effects that break the detector sym- 
metry. We therefore define v as the average of the values 
of V from the sixteen proton-cells at the same |z|, i.e. 
z = ±2, ±6, ±10, and ±14 cm. Each set of detectors 
corresponds to a paired-ring with the same symmetry as 
the full detector, e.g. the shaded detectors in Fig. [TJ We 
then define 



where i^z? = 0.378 is the average of f- (Kg''" -Kg de- 
termined by Monte Carlo. The experiment provides four 
independent measurements corresponding to each of the 
four paired-rings. The systematic corrections to D pre- 
sented in Table I yield our final value for D. Eqn. [S]is 
based on the following: 1) accurate background correc- 
tions, 2) uniform proton and electron detection efficien- 
cies, 3) cylindrical symmetry of the neutron beam and 
polarization, and 4) accurate determination of Kb^ P, 
and spin state. 

Backgrounds not properly accounted for contribute 
two systematic errors: 1) multiplicative errors due to di- 
lution of the asymmetries, and 2) spin-dependent back- 
grounds that can lead to a false D. Errors in background 
subtraction, as well as possible spin-dependent asymme- 
tries in this background, have a small effect. The multi- 
plicative correction to the value of w^^^^ due to backscat- 
tered electrons was determined using the Monte Carlo. 



^ Zero indicates no correction applied. ^ Includes spin-flip 
time, cycle asymmetry, and flux variation. Included in the 
definition of D. Assumed polarization uncertainty of 0.05. 

The uncertainty given in Table |T] refiects the 20% un- 
certainty assigned to the backscattering fractions due to 
limitations of the detector and beam model and due to 
limited knowledge of backscattering at energies below a 
few hundred keV. Proton backscattering, though observ- 
able, produces a negligible effect on D. 

In principle, the values of w^''^^ are independent of the 
absolute efficiencies of the proton and electron detectors; 
however, they do depend on any energy dependence of 
the efficiencies through the factors (pe/Ee). Spatial vari- 
ation of the efficiencies breaks the symmetry assumed in 
combining proton-cell data into paired-rings. Beta en- 
ergy thresholds were observed to vary less than 20 keV 
across the detector, implying the almost negligible cor- 
rection given in Table HI Proton detector efficiency varia- 
tions, however, were more significant. Lower energy pro- 
ton thresholds varied across the detector and over the 
course of the experiment. These thresholds combined 
with the spin-dependence of the accelerated proton en- 
ergy spectra can result in significant deviations in the 
values of w^''^^ , though the effect on the the value of 
is largely mitigated because the low-energy portion 
of the proton energy spectrum is roughly the same for 
the cr and cl coincidence pairs. To estimate the proton- 
threshold- nonuniformity effect on D, spin-dependent pro- 
ton energy spectra were generated by Monte Carlo for 
all proton-detector-electron-detector pairings and convo- 
luted with model detector response functions based on 
fits to the average proton-SBD spectra. The average 
fit parameters were varied over a range characteristic of 
the observed variations during the run. Representative 
thresholds were then applied to determine the effect on 
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D. An alternative and consistent estimate was derived by 
correcting the values of w^'^'^^ on a day- by-day basis using 
the spectrum centroid shift and empirically determined 
functional form of the spectrum at the threshold. 

Beam expansion from a radius of 2.5 cm to 2.75 cm 
combined with the magnetic field breaks the symmetry of 
the detector because the average proton-electron opening 
angle for each proton-electron detector pair is modified. 
Monte Carlo calculations using measured upstream and 
downstream density profile maps were used to calculate 
the correction given in Table HI Possible inaccuracies in 
the determination of the beam density were estimated 
and their implications explored with the Monte Carlo 
code. The effect on the value of v^^ as a function of 
position is illustrated in Fig. |3l 
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FIG. 3: Solid (open) squares show the values of v averaged 
over the four planes for proton cells on the even (odd) side of 
the proton detection plane (the side with P2 vs Pi in Fig.[l|. 
Monte Carlo results are indicated by lines. The broken sym- 
metry, due to the combination of the magnetic field and beam 
expansion, is evident in the shift in the crossing point from 
the detector center. 



spin-flipper, as well as the uncertainty in the instrumen- 
tal constant Ke,, can be shown to produce asymmetries 
proportional to D. These effects are listed in Table HI 

Correlations of D with a variety of experimental pa- 
rameters were studied by varying the cuts and by break- 
ing the data up into subsets taken under different con- 
ditions of proton acceleration voltage and number of live 
SBDs as shown in Fig. |4l A linear correlation of D with 
high-voltage, revealed by the cuts study, yields — 5.6 
with 11 DOF compared to 10.4 for 12 DOF for no cor- 
relation. The acceleration-voltage dependence of the fo- 
cusing properties was extensively studied by Monte Carlo 
with no expected effect, and we intrepret the 2.1 sigma 
slope as a statistical fluctuation. 

To improve the symmetry of the the detector, we com- 
bine counts for the entire run to determine the values of 
and to extract D for each paired-ring. A blind analy- 
sis was performed by adding a constant hidden factor to 
Eqn. [21 The blind was only removed once all analyses of 
systematic effects were complete and had been combined 
into a final result. The weighted average of D for the four 
paired-rings is 0.72 ± 1.89 with = 0-73 for 3 DOF. 
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For a symmetric beam, contributions to v^^ due to 
transverse polarization cancel for opposing proton planes; 
however small azimuthal beam asymmetries can affect 
this cancellation. This asymmetric-beam transverse- 
polarization (ATP) effect is proportional to both sin 6'p, 
the angle of the average neutron-spin orientation with 
respect to the detector axis, and sin(0p — ^''•), where 
(ff^ is the effective azimuthal position of the proton cell, 
and (pp \s the azimuthal direction of the neutron polar- 
ization. To study this effect transverse-polarization cal- 
ibration runs with 9p = 90° and several values oi 4>p 
were taken over the course of the experiment. In these 
runs the ATP effect was amplified by « 200. The val- 
ues of sin^p and (j)p for the experiment were determined 
using the calibration runs and Monte Carlo corrections 
for the beam density variations. To estimate the effect, 
the extreme value of suvOp — 12.8 x 10^'^ and the range 
of —31.5° < (pp < 112.2° were used. The uncertainty is 
due to uncertainties in the angles 9p and (pp. The effect 
of nonuniform beam polarization is also given in Table ID 
Time-dependent variations in flux, polarization, and the 



FIG. 4: Results for D by run subset. Uncertainties are statis- 
tical. The fully functioning paired-rings used for each subset 
are indicated: 1-4 indicates all paired-rings were used. 

The result including all corrections from Table U is 

D = (-0.96 ± 1.89{stat) ± 1.01(s?;s)) x lO"'^. 

Our result represents the most sensitive measurement of 
the D coefficient in nuclear beta decay. Assuming purely 
vector and axial-vector currents, (/)yiy = 180.013°±0.028° 
which is the best direct determination of a possible CP- 
violating phase between the axial and vector currents in 
nuclear beta decay. Previously the most sensitive mea- 
i9Ne, 
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surement was in ^"Ne, with D = (1 ± 6) x lO"'* 
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